Algebraic methods for counting Euclidean 
embeddings of rigid graphs 



Ioannis Z. Emiris 1 , Elias P. Tsigaridas 2 , and Antonios Varvitsiotis 

1 National and Kapodistrian University of Athens, Greece. 
2 INRIA Mediterranee, Sophia- Antipolis, France. 



Abstract. The study of (minimally) rigid graphs is motivated by numerous 
applications, mostly in robotics and bioinformatics. A major open problem con- 
cerns the number of embeddings of such graphs, up to rigid motions, in Eu- 
clidean space. We capture embeddability by polynomial systems with suitable 
structure, so that their mixed volume, which bounds the number of common 
roots, to yield interesting upper bounds on the number of embeddings. We 
focus on R 2 and R 3 , where Laman graphs and 1-skeleta of convex simplicial 
polyhedra, respectively, admit inductive Henneberg constructions. We establish 
the first lower bound in R 3 of about 2.52", where n denotes the number of ver- 
tices. Moreover, our implementation yields upper bounds for n < 10 in R 2 and 
R 3 , which reduce the existing gaps, and tight bounds up to n = 7 in R 3 . 

Keywords: rigid graph, Euclidean embedding, Henneberg construction, poly- 
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1 Introduction 

Rigid graphs (or frameworks) constitute an old but still very active area of research due 
to their deep mathematical and algorithmic questions, as well as numerous applications, 
e.g. mechanism and linkage theory [12, 13], and structural bioinformatics [6, 9, 11]. 

Given a graph G = (V, E) and a collection of edge lengths kj £ R+, for e E, 
a Euclidean embedding in R d is a mapping of V to a set of points in R d , such that 
lij equals the Euclidean distance between the images of the i-th and j-th vertices, for 
€ E. Euclidean embeddings impose no requirements on whether the edges cross 
each other or not. A graph is (generically) rigid in R d iff, for generic (or random) edge 
lengths, it is embedded in R d in a finite number of ways, modulo rigid motions. A graph 
is minimally rigid iff it is no longer rigid once any edge is removed. 

A graph is called Laman iff \E\ ~ 2\V\ — 3 and, additionally, all of its induced 
subgraphs on k < \V\ vertices have < 2k — 3 edges. The class of Laman graphs co- 
incides with the generically minimally rigid graphs in R 2 , and also admit inductive 
constructions. In R 3 there is no analogous combinatorial characterization of generically 
rigid graphs. On the other hand, the 1-skeleta, or edge graphs, of (convex) simplicial 
polyhedra are minimally rigid in R 3 , and admit inductive constructions, cf. Section 4. 

In this paper, we deal with the problem of computing the maximum number of 
distinct planar and spatial Euclidean embeddings of (minimally) rigid graphs, up to 
rigid motions, as a function of the number of vertices. To study upper bounds, we 
define a square polynomial system, expressing the edge length constraints, whose real 
solutions correspond precisely to the different embeddings. Here is a system expressing 
embeddability in R 3 , where (xi,yi,Zi) are the coordinates of the i-th vertex, and 3 



vertices are fixed to discard translations and rotations: 

f x t = at, yt = h, Zi = c», i = l,2,3, a^&i^el, , . 

\(^-^) 2 + (y l -y J ) 2 + (z l -^-) 2 = ;? J , (i,j) {(1,2), (1,3), (2,3)} 1 j 

All nontrivial equations are quadratic; there are In — 4 for Laman graphs, and in — 9 
for 1-skclcta of simplicial polyhedra, where n is the number of vertices. The classical 
Bezout bound on the number of roots equals the product of the polynomials' degrees, 
and yields 4™~ 2 and 8™~ 3 , respectively. It is indicative of the hardness of the problem 
that efforts to substantially improve these bounds have failed. 

For the planar and spatial case, the best known upper bounds are ( 2 ™JT 2 4 ) ~ 4"~ 2 /y / 7r(n 

and iB^^n-z) ~ 8™~ 3 /((n - 2) v / 7 r(n- 3)), respectively. These bounds were ob- 
tained using complex algebraic geometry [1,2]. For R 2 , there exist lower bounds of 
24 L("-2)/4J _ 2 _2i« an d 2 • l2L(«-3)/3j ~ 2.29™/6, obtained by a caterpillar and a 
fan 3 construction, respectively. Both of these bounds are based on the Dcsargues (or 
3-prism) graph (Figure 2). 

In applications, it is crucial to know the number of cmbcddings for specific (small) 
values of n. The most important result in this direction was to show that the Desargues 
graph admits 24 embeddings in the plane [2]. Moreover, the K33 graph admits 16 
embeddings in the plane [14, 12] and the cyclohexane graph admits 16 embeddings in 
space [6]. 

Mixed volume (or Bernstein's bound) of a square polynomial system exploits the 
sparseness of the equations to bound the number of common roots, it is always bounded 
by Bezout's bound and typically much tighter, cf. Section 3. We have implemented 
specialized software that constructs all rigid graphs up to isomorphism, for small n, 
and computes the mixed volumes of their respective polynomial systems. Our results 
indicate that mixed volume can be of general interest in enumeration problems. 

Our main contribution is twofold, besides some straightforward upper bounds in 
Lemmas 2 and 7. First, we derive the first lower bound in R 3 : 

16 L(«-3)/3j -2.52", n> 9, 

by designing a cyclohexane caterpillar. 

Second, we obtain upper and lower bounds for n < 10 in R 2 and R 3 , which reduce 
the existing gaps, see Tables 1 and 2 in the Appendix. Moreover, we establish tight 
bounds up to n = 7 in R 3 by appropriately formulating the polynomial system. We 
apply Bernstein's Second theorem to show that the naive polynomial system cannot 
yield tight mixed volumes, a fact already observed in the planar case [10]. 

The rest of the paper is structured as follows: Section 2 discusses the case d = 2, 
Section 3 presents our algebraic tools and our implementation, Section 4 deals with R 3 , 
and we conclude with open questions and a conjecture. Omitted Figures and Tables 
have been included in the Appendix. 

Some results appeared in [7] in preliminary form. 

2 Planar embeddings of Laman graphs 

Laman graphs admit inductive constructions that begin with a triangle, followed by 
a sequence of Henneberg-1 (or Hi) and Hennenerg-2 steps (or H 2 ). Each such step 
adds a new vertex as follows: a Hi step connects it to two existing vertices; a H 2 



3 This corrects the exponent of the original statement. 



step connects it to three existing vertices having at least one edge among them, which 
is removed (Figure 1). We represent each Laman graph by As4...,s„, where Si € 
{1,2}; this is known as its Henneberg sequence. A Laman graph is called Hi iff it 
can be constructed using only Hi steps, and H 2 otherwise. Since two circles intersect 
generically in two points, a Hi step at most doubles the number of embeddings and this 
is tight, generically. It follows that a Hi graph on n vertices has 2™~ 2 embeddings [2]. 

One can easily verify that every A 2 graph is isomorphic to a Al graph and that 
every A 12 graph is isomorphic to a All graph. Consequently, all Laman graphs with 
n = 4, 5 are Hi and they have 4 and 8 embeddings, respectively. 

For a Laman graph on 6 vertices, there are 3 possibilities: it is either Hi, or 
the Desargues graph. Since the ^3,3 graph has at most 16 embeddings [14, 12] and the 
Desargues graph has 24 embeddings [2], the uppper bound is 24 for n = 6. 

Lemma 1. The maximum number of Euclidean embeddings for Laman graphs with 
n = 7, . . . , 10 is 64, 128, 512 and 2048, respectively. 

Proof. Using our software (Section 3), we construct all Laman graphs with n = 7, . . . , 10, 
and compute their respective mixed volumes, thus obtaining the upper bounds. 

Table 1 in the Appendix summarizes our results for n < 10. The lower bound for 
n = 9 follows from the Desargues fan [2]. All other lower bounds follow from the fact 
that a Hi step exactly doubles the number of embeddings. 

We now establish an upper bound, which improves upon the existing ones when 
our graph contains many degree-2 vertices. Our proof parallels that of [10] which uses 
mixed volumes to bound the effect of a Hi step, when it is the last one in the Henneberg 
sequence. 

Lemma 2. Let G be a Laman graph with k > 4 degree-2 vertices. Then, the number 
of planar embeddings of G is bounded above by 2 k ~ 4 4 n ~ k . 

Proof. The removal of the existing k degree-2 vertices does not affect any of the ex- 
isting degree-2 vertices (because the remaining graph is also Laman), although it may 
create new ones. Since the remaining graph has n — k vertices, the Bezout bound of 
its polynomial system is equal to 4 n ~ k and thus the number of embeddings is bounded 
above by 2 fc - 4 4™-' £ . 

3 An algebraic interlude 

This section introduces mixed volumes and discusses our computer-assisted proofs; for 
background see [4, 5] and references therein. 

Given a polynomial / in n variables, its support is the set of exponents in N™ 
corresponding to nonzero terms (or monomials). The Newton polytope of / is the 
convex hull of its support and lies in M. n . Consider poly topes Pi C W 1 and parameters 
Aj e K, Aj > 0, for i = 1, . . . , n. Consider the Minkowski sum of the scaled polytopes 
A1P1 + • • • + \ n P n € K"; its (Euclidean) volume is a homogeneous polynomial of 
degree n in the Ai. The coefficient of Ai • • • A„ is the mixed volume of Pi, . . . ,P n . If 
Pi = ■ ■ ■ = P n , then the mixed volume is nl times the volume of Pi. We focus on the 
topological torus C* = C - {0}. 

Theorem 3. [4] Let fi = •■■ = /« = be a polynomial system in n variables with real 
coefficients, where the fi have fixed supports. The number of isolated common solutions 
in (C*) n is bounded above by the mixed volume of (the Newton polytopes of) the fi. 
This bound is tight for a generic choice of coefficients of the fi 's. 



Bernstein's Second Theorem below, describes genericity. Given v € R" — {0} and 
polynomial /,, we denote by d v fi the polynomial obtained by keeping only those terms 
whose exponents minimize the inner product with v. The Newton polytope of d v fi is 
the face of the Newton polytope of fi supported by v. 

Theorem 4. [4] If for all v e K" — {0} the face system d v fi = . . . = d v f n = has 
no solutions in (C*)™, then the mixed volume of the fi exactly equals the number of 
solutions in (C*)™, and all solutions are isolated. Otherwise, the mixed volume is a 
strict upper bound on the number of isolated solutions. 

This theorem was used to study planar embcddings [10]; we shall apply it to R 3 . 

Now, we describe how we make use of mixed volume in our implementations. In 
order to bound the number of embeddings of rigid graphs, we have developed specialized 
software that constructs all Laman graphs and all 1-skeleta of simplicial polyhedra in 
R 3 with n < 10. Our computational platform is SAGE 4 . We exploit the fact that these 
graphs admit inductive constructions, and construct all graphs using the Henneberg 
steps. The latter were implemented, using SAGE's interpreter, in Python. After we 
construct all the graphs, we classify them up to isomorhism using SAGE's interface 
for N.I.C.E., an open-source isomorphism check engine, keeping for each graph the 
Henneberg sequence with largest number of Hi steps. 

Then, for each graph we construct a polynomial system whose real solutions express 
all possible embeddings, using formulation (2) below. For each system we bound the 
number of its (complex) solutions by computing its mixed volume, using [5]. Notice 
that, by genericity, solutions have no zero coordinates. For every Laman graph, to 
discard translations and rotations, we assume that one edge is of unit length, aligned 
with an axis, with one of its vertices at the origin. In R 3 , a third vertex is also fixed 
so as to belong to a coordinate plane. The corresponding coordinates are given specific 
values and are no longer unknowns. Depending on the choice of the fixed edge, we 
obtain different systems hence different mixed volumes. Since they all bound the actual 
number of embeddings, we use the minimum mixed volume. 

We used an Intel Core2, at 2.4GHz, with 2GB of RAM. We tested more that 20, 000 
graphs and computed the mixed volume of more than 40, 000 polynomial systems. The 
total time of experiments was about 2 days. Tables 1 and 2 summarize our results. 

4 Spatial embeddings of 1-skeleta of simplicial polyhedra 

This section extends the previous results to 1-skeleta of (convex) simplicial polyhedra, 
which are minimally rigid in R 3 [8]. For such a graph (V,E), we have \E\ = 3\V\ — 6 
and all of the induced subgraphs on k < \V\ vertices have < 3k — 6 edges. 

Consider any k + 2 vertices forming a cycle with > k — 1 diagonals, k > 1. The 
extended Henneberg-fc step (or Hk), k = 1,2,3, corresponds to adding a vertex, con- 
necting it to the k + 2 vertices, and removing k — 1 diagonals among them, cf. Figure 4. 
It is known that a graph is the 1-skeleton of a simplicial polyhedron in R 3 iff it has 
a construction that begins with the 3-simplex followed by any sequence of Hi,H 2 , H 3 
steps [3]. 

Since 3 spheres intersect generically in two points, a H± step at most doubles the 
number of spatial embeddings and this is tight, generically. In order to discard trans- 
lations and rotations, we fix a (triangular) facet of the polytope; we choose wlog the 



4 http : / /www . sagemath . org/ 



first 3 vertices and obtain system (1) of dimension 3n. It turns out that this system 
does not capture the structure of the problem. 

Specifically choosing direction v = (0, 0, 0, 0, 0, 0, 0, 0, 0, -1, . . . , -1) e R 3 ™, the 
corresponding face system is: 



This system has (ai, b\, ci, . . . , 03, 63, C3, 1, 1, 7\/2, ...,1,1, 7V2) <= (C*) 3 ™ as a solution, 
where 7 = ±-\/— 1. According to Theorem 4, the mixed volume is not a tight bound on 
the number of solutions in (C*) 3 ". This was also observed, for the planar case, in [10]. To 
remove spurious solutions (at toric infinity), we introduce variables Wi — xf + yf + zf, 
for i — l,...,n. This yields the following equivalent system, but with lower mixed 
volume: 



This is the formulation we will use for our computations. 

For n = 4, the only simplicial polytope is the 3-simplex, which admits 2 embeddings. 
For n = 5, there is a unique graph that corresponds to a 1-skeleton of a simplicial 
polyhedron, see Figure 6 [3]. This graph is obtained from the 3-simplex trough a Hi 
step, so for n = 4 the bound is 4 and it is tight. 

Theorem 5. The 1-skeleton of a simplicial polyhedron on 6 vertices has at most 16 
embeddings and this bound is tight. 

Proof. There are two non-isomorphic graphs Gi, G2 for n — 6 [3], see Figure 3. For G\, 
the mixed volume is 8. Since all facets of G2 are symmetric, we fix one and compute 
the mixed volume, which equals 16, so the upper bound is 16. G2 is the graph of the 
cyclohexane, which admits 16 different Euclidean embeddings [6]. To see the equiva- 
lence, recall that the cyclohexane is essentially a 6-cycle (fig. 5) with known lengths 
between vertices at distance 1 (adjacent) and 2. The former are bond lengths whereas 
the latter are specified by the angle between consecutive bonds. 

We now pass to general n and establish the first lower bound in M. 3 . 

Theorem 6. There exist edge lengths for which the cyclohexane caterpillar construc- 
tion has i6LC«-3)/3j _ 2.52" embeddings, for n > 9. 

Proof. We glue together copies of cyclohexanes sharing a common triangle. The result- 
ing graph is the 1-skeleton of a simplicial polytope. Each new copy adds 3 vertices, and 
since there exist edge lengths for which the cyclohexane graph has 16 embeddings the 
claim follows. An example with 2 copies is illustrated in Figure 7. 

Table 2 summarizes our results for n < 10. The upper bounds for n = 7, . . . , 10 are 
computed by our software. The lower bound for n = 9 follows from the cyclohexane 
caterpillar. All other lower bounds are obtained by applying a Hi step to a graph with 
one fewer vertex. Lastly, we state without proof a result similar to Lemma 2. 

Lemma 7. Let G be the 1-skeleton of a simplicial polyhedron with k > 9 degree-3 
vertices. Then the number of embeddings of G is bounded above by 2 k ~ 9 8 n ~ k . 





(2) 



5 Further work 



Undoubtedly, the most important and oldest problem in rigidity theory is the com- 
binatorial characterization of rigid graphs in R 3 . In the planar case, existing bounds 
are not tight despite our earlier attempt to close all gaps up to n — 8. This is due 
to the fact that root counts include rotated copies of certain embeddings. We expect 
that a more careful approach should settle these cases. Since we deal with Henneberg 
constructions, it is important to determine the effect of each step on the number of 
embeddings: a Hi step always doubles their number; we conjecture that H 2 multiplies 
it by < 4 and spatial H3 by < 8, but these may not always be tight. Our conjecture 
has been verified for small n. Here, mixed volume may help; it is also challenging to 
understand its relevance, despite the fact it is a bound on complex roots. As for lower 
bounds, specific graphs with a large number of embeddings may be combined to yield 
tighter bounds. 
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Table 1. Bounds and Henneberg sequences for Laman graphs for n < 10. Bold text indicates 
the Henneberg sequence of the graph yielding the upper bound. 



n = 


4 


5 


6 


7 


8 


9 


10 


lower 


2 


4 


16 


32 


64 


256 


512 


upper 


2 


4 


16 


32 


160 


640 


2560 


Hi 


A 


Al 


All 


Alll 


Al 4 


Al° 


Al b 








A12 


Al 2 2 


Al 3 2 
Al 2 2 2 


Al 4 2 
Al 3 2 2 


Al & 2 
Al 4 2 2 


H 2 










Al 2 21 


Al 3 21 
Al 2 2 3 


Al 4 21 

Al 3 2 3 

Al 3 21 2 

Al 3 212 

Al 3 2 2 l 

Al 2 2 4 


H 3 

















Table 2. Bounds and Henneberg sequences for 1-skeleta of simplicial polyhedra for n < 10, 
where A is the 3-simplex. H3 need not apply before n — 13. Bold text indicates the Henneberg 
sequence of the graph that yields the upper bound. 
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Fig. 1. The planar Henneberg steps; the bot- 
tom vertex is new. 
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Fig. 2. The Desargues graph. 




Fig. 6. The only 1-skeleton of a simplicial poly- Fig. 7. A Cyclohexane caterpillar with 2 
tope on 5 vertices. copies. 



